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I. INTRODUCTION 

Almost a century after it was first proposed mm and fifty years since the invention of the laser allowed its first 
observation [3], the phenomenon of stimulated Brillouin scattering (SBS) may only now be entering its golden age. 
At its simplest, SBS refers to the stimulated interaction between a pair of coherent optical waves and a resonant 
hypersonic acoustic wave. SBS has traditionally been encountered as the scattering of an optical pump beam into 
a backward traveling Stokes beam of slightly lower frequency by an acoustic wave oscillating at the optical beat 
frequency. The acoustic wave is generated by the process of electrostriction 015], and both the Stokes and acoustic 
wave grow by a process of positive feedback. In optical fiber, this process can be highly efficient and is often 
described as the “strongest” fiber nonlinearity. This can be problematic, as SBS prevents the propagation of high 
power narrow-bandwidth pumps. Nevertheless, SBS in fibers has long provided a mechanism for producing narrow 
linewidth lasers and amplifiers IMS], filters and other spectral components for microwave photonics [sini, as well as 
various sensors min]. 

Like most optical nonlinearities, the emergence of sub-wavelength scale waveguides with strong confinement and 
tunable dispersion has greatly increased the efficiency, utility and reach of stimulated Brillouin processes. This 
includes the generation of frequency combs by cascaded Brillouin generation in microstructured small-core fibers in 
both backward m and forward [16] configurations, slow and fast light effects and novel microstructured hber 
lasers [lids]. 

Motivated by such studies, the development of on-chip SBS in highly-nonlinear integrated waveguides has recently 
been pursued aggressively |5D]. Attaining efficient on-chip SBS is complicated by the requirement of simultaneous 
confinement of both the optical and acoustic fields. This is non-trivial because optically dense materials suitable 
for optical waveguide cores are commonly mechanically stiff and therefore susceptible to leakage of the acoustic 
wave into softer substrates; the silicon on silica system is an important example. Consequently, on-chip SBS was 
first achieved m in rib waveguides made from nonlinear chalcogenide glasses, which combine high refractive index 
and nonlinearity with relative mechanical softness. Subsequently, SBS in silicon waveguides has been observed in 
Si/SiN membranes [H] and elevated rails [531 dl], which both exploit physical isolation of the waveguide to minimize 
acoustic losses. Considerable development will be needed to reach designs suitable for mass-fabrication, but a practical 
platform for on-chip SBS would enable numerous applications |25j in microwave photonics nni [niiMii], sensing, 
isolators |30l |3T] and chip-based lasers |32l [33] . 

A key driver for developing SBS in sub-micron waveguides was the realization by Rakich et al. |34l [35] that 
at small scales, there are new contributions to SBS associated with radiation pressure of light on the waveguide 
boundaries, and the back-action of “moving boundaries” on the optical field [33]. Depending on the particular 
waveguide configuration and combination of optical modes, these contributions can either reinforce or counteract 
the more familiar bulk contributions from electrostriction and photoelasticity [34]. As well, waveguides in which the 
acoustic fields are strongly confined can enhance the scattering efficiency of near-stationary quasi-transverse acoustic 
waves, a requirement for efficient forward SBS where the pump and Stokes wave co-propagate [MISS]. 

Following this realization there was some variation in the literature as to how best to incorporate the new effects into 
a coupled mode theory self-consistently. In conventional SBS, electrostriction (the mechanical stress induced by the 
optical field) is accompanied by the complementary process of photoelasticity (the change in the dielectric response 
induced by the acoustic strain) [3]. The two processes are captured by identical coupling terms in the coupled mode 
theory, as required by the Manley-Rowe relations [531137]- One should expect the same symmetry between the 
effect of radiation pressure on the waveguide boundaries driving the acoustic field, and the reverse effects of moving 
boundaries on the optical field. However, the initial formulations in terms of optical forces and the Maxwell stress 
tensor led to some confusion about whether certain additional coupling terms arise or whether they are essentially 
“double-counting”. These include concepts of an electrostrictive “boundary pressure” [55] and bulk contributions to 
the radiation pressure term [311 ES]. 

Recently, our group provided a new derivation of the coupled mode equations m that avoids the formalism of optical 
forces. Instead we used thermodynamic arguments to unambiguously identify the correct coupling term describing 
both radiation pressure and moving boundaries effects. We found that this term indeed matches expressions suggested 
in several papers [341138j and the matter seems to be resolved. It turns out, for instance, that the appearance of 
an electrostrictive boundary pressure depends on whether electrostriction is viewed in terms of a stress or a force 
density [37]. Nevertheless, the argument establishing the correct form for the moving boundary coupling was quite 
involved, and a simpler, more direct derivation would be desirable. 

In this work, we provide such a derivation. Rather than the standard approach of applying slowly-varying envelope 
approximations to the wave equation, we extend a quantized multimode Hamiltonian formalism of integrated optical 
waveguides [39] to include opto-acoustic interactions. The photoelastic and radiation pressure couplings are introduced 
through a single interaction energy term in the Hamiltonian, and fully quantum equations of motion are obtained 
from the Heisenberg equations. In the classical limit, the coupled mode equations of Wolff et al. m emerge naturally 
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and simply, with no ambiguity about double-counting. 

There are number of other advantages to our approach. The fundamental quantum process underlying SBS—the 
stimulated decay of a pump photon into a lower energy Stokes photon and an acoustic phonon—is manifestly visible 
in the interaction term. Further, obtaining quantum equations that respect the appropriate operator commutation 
relations provides an important starting point for the investigation of effects at the boundary of quantum and classical 
opto-acoustics. This is likely to become more important as the distinction between cavity optomechanics and guided 
wave opto-acoustics becomes increasingly blurred snun], and phonon confinement strategies are improved. The 
propagation of guided wave acoustic fields with strongly modified phonon density of states is likely not far off, which 
raises the prospect of Brillouin interactions in the quantum regime. Finally, the Hamiltonian formalism has proved 
very powerful for the description of other integrated quantum nonlinear processes such as spontaneous four wave 
mixing |42L I43j and spontaneous parametric downconversion [44] . 

We should note that quantum or analytical dynamics approaches to guided wave opto-acoustics themselves have 
some pedigree. Hamiltonian approaches to SBS date back to the first rigorous treatment by Shen and Bloembergen 
who gave an analysis for plane waves in the semi-classical limit. Drummond and Corney incorporate Raman gain 
into their quantized theory of nonlinear fiber propagation |45j . In that case, the Raman response, which depends on 
the detailed glass composition and network, is introduced through a phenomenological measured response function. 
In contrast, for the Brillouin couplings we consider here, the phonon response is entirely determined by the bulk 
elastic properties and waveguide geometry, and so can be calculated exactly using acoustic mode solvers. Finally, van 
Laer et al. |40| have recently discussed the connections between quantum optomechanics for single or few resonator 
systems and classical SBS [ID] . They identify an elegant connection between the opto-acoustic coupling in waveguide 
SBS and the corresponding coupling in quantum optomechanical systems. The latter is treated with a single mode 
Hamiltonian approach which is appropriate for the optomechanics of a resonator consisting of a single cavity, but 
limits its application to longer structures with continuous phonon spectra. In contrast, for the waveguide problem 
that is our focus, a full multi-mode treatment is appropriate in order to treat arbitrary input optical fields. 

The paper is structured as follows. In Section [ll| we construct a Hamiltonian description of acoustics, including 
useful expressions for group velocity and power flow. In Section m we review some necessary results from guided 
wave electromagnetic quantization. In Section EYl which is the core of the paper, we provide the full optoacoustic 
Hamiltonian, find expressions for the coupling terms, and show directly how the symmetry between radiation pressure 
and moving boundary effects emerges. In Section]^ we derive quantum coupled mode equations for the system, make 
connections to prior expressions for the coupling strength, and recover the classical coupled mode equations for SBS. 
Finally, in Section jVTj we discuss directions for future examination including the important issue of phonon dissipation. 
A comprehensive supplementary materials document provides detailed derivations of many of the results. 


II. HAMILTONIAN FORMULATION OF GUIDED WAVE ACOUSTICS 

The classical theory of guided elastic waves is of course very mature and can be formulated in many guises. Auld [15] 
provides an excellent introduction for readers with an optics background. Since our goal is a Hamiltonian operator 
representing the complete opto-acoustic system we begin by re-framing guided acoustic wave propagation in a quantum 
Hamiltonian picture; we have not found such a formulation in the literature. 

For purely classical applications, one might build a Hamiltonian from which the dynamics are determined by 
Hamilton’s equations. For generality, we construct the theory in a quantum form, using canonical quantization with 
commutators that follow from the standard association with the Poisson brackets of the classical formulation: 


A. Hamiltonian operator 

To identify a classical theory suitable for canonical quantization we should begin with canonical variables, which 
will become the canonical operators in the quantum theory, and a classical Hamiltonian that both yields the standard 
equations of motion in the form of the elastic wave equation, and is numerically equal to the classical energy of the 
system. To that end we introduce vector held variables u(r) describing the displacement and 7r(r) as their conjugate 
momenta; these will become operators in the quantum theory, although we will not explicitly include “hats” in our 
notation. 

Following standard quantum mechanics, we naturally choose the commutation relations 


[zi”(r), 7r'"(r')] = - r'). 


( 2 ) 
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FIG. 1: Schematic geometry of the interacting waves in integrated SBS. Arrows at the top denote wavenumbers for 
forward SBS; arrows at the bottom denote wavenumbers for backward SBS. 


where superscripts stand for Cartesian components. The obvious acoustic Hamiltonian operator is 

H^ = J ^ J S^Hry^>^\r)S'^\r)dr, (3) 

where the integration is over all space. Here p{r) and are position-dependent c-number quantities representing 

the density and stiffness tensor respectively, while 




(4) 


is the strain tensor operator [45], and again we emphasise that unless stated otherwise, u(r) and 7r(r) are to be read 
as operators. In (|^ and throughout, repeated indices are to be summed over the Cartesian coordinates x, y, z. We 
will formally assume these quantities are continuous functions of position, although they may change drastically as 
one moves from solid to air, for example. The limit to these functions changing discontinuously can be taken at the 
end of the calculation when matrix elements and the like are evaluated, but in all dynamical equations and derivations 
p{r) and should be taken as continuous functions. 

Let us comment on the precise meaning of the operators u and tt. In elastic theory, the displacement field applies 
to “volume elements” much larger than the atomic scale, but much smaller than the characteristic wavelength of 
any acoustic excitation. Therefore, the displacement and momentum operators do not correspond directly to any 
individual physical oscillator but describe collective excitations of the mesoscopic bulk medium. Nevertheless the 
low-energy phonon excitations that emerge from the theory are very real and their quantization is physical. For 
example, completeness relations are correct up to an appropriate wavevector cutoff, far above the typical wavevector 
of any Brillouin-induced excitation. 

The strain tensor Q is obviously symmetric in its two indices, and since the stiffness tensor appears with two strain 
tensors in ([^, c®^^^(^can be taken to satisfy as well as In all then [47] , 


so we can write ([^ as simply 


7r*(r)7r*(r) 

2p(r) 



dd dr'^ 


dr. 


Evaluating the Heisenberg equations of motion (see supplementary material section S.I) 


= [7r"(r),77^] , 


( 6 ) 


(7a) 

(7b) 
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yields 




m 


(r.O 


7r"(r,i) 


p(r) 


(8a) 

(8b) 


Classically these equations indicate, as expected, that the conjugate momentum is the product of the mass density 
and the velocity field u", and that the momentum evolves according to Newton’s laws. 

From Q, we recover the standard acoustic wave equation [35]: 


P{^) 


i9^M”(r,t) 


J)_ 

dri 



(9) 


Thus we hav e th e desired equations of motion in both the classical and quantum descriptions. Moreover, substituting 
the relation (8a I into ^ and taking the classical limit by dropping the operator nature of u and , we recover the 
classical Hamiltonian 




T + V. 


This is clearly numerically equal to the sum of kinetic and potential energy |46j . and so @ meets all the criteria 
to be an appropriate Hamiltonian for quantization. Note that at this point, no acoustic dissipation (viscosity) has 
been included. Direct incorporation of losses is less straightforward in Hamiltonian approaches than in Lagrangian 
approaches, in which a Rayleigh dissipation function can be introduced. The phonon loss is important to the physics 
of SBS, but can be incorporated later perturbatively. We discuss this in Section lYg 


B. Modes and new fields 


We now seek to reduce the Hamiltonian to standard harmonic oscillator form; the low energy states will describe 
the quasi-particle or phonon excitations of the medium. To diagonalise the acoustic Hamiltonian we look for classical 
solutions of 0 of the form 

u(r,t) =UA(r)e-*^^*+c.c., (10) 

where A is a mode index. This implies 

- ^ ( 11 ) 


where again we have used the fact that c”'’^^(r) is symmetric in its last two indices. It 
linear operator generating mode functions is Hermitian, which the form in Eq. (Ill is 
system, we introduce new fields 


is convenient however, if the 
not. To obtain a Hermitian 


u(r) = y/p{r) u(r) 

*(r) = 




which clearly preserve the commutation relations. That is, 

[w”(r),if'"(r')] =iM”™5(r-r'). 

We now look for mode solutions in terms of the new fields in the form 

u(r,t) =UA(r)e-*^^‘ + c.c., 
ff(r,t) =nA(r)e-*^''*+c.c. 


( 12 ) 


(13) 


(14a) 

(14b) 
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We introduce the operator which acts on a general vector function C(r) as 

1 d .d(C^ 


7W”'=(r)C"'(r) = - 


yp(r)5rJ 







(15) 


Section S.II of the supplementary material shows that A^"^(r) is indeed Hermitian. Then, using integration by parts 
we can write the Hamiltonian (|^ in terms of the new fields (12) in the simple form 


= o / 7r®(r)7r*(r) dr + - / it*(r)Al**^(r)u*^(r) dr. 


(16) 


As the basis of our acoustic modal expansion, Eq. (151 plays an analogous role to that which the vector Helmholtz 
equation (commonly known as the “master equation” in the photonic crystal literature [48],) plays in the quantization 
of the electromagnetic field (see (43) below). Its Hermitian form is useful both for developing the formalism but also 
for formulating mode-solving algorithms based on energy functional minimization |48) . 

Now if we construct eigenfunctions FA(r) of such that 




it follows from 0 and ([T^ that we can take 


UA(r) = FA(r), 
f[A(r) = -zHAFA(r), 


(17) 

(18) 


where the second equation follows from substituting (14) in (8a). 

From the Hermiticity of Al”^, the eigenfunctions FaJt) with different eigenvalues will automatically be orthogonal, 
and if there are eigenfunctions with the same eigenvalue we can construct them to be orthogonal. Thus we can write 


J Fl(r).FA'(r) dr = (5AA'. 


(19) 


Here the integration is over a finite volume that, in the end, can be allowed to pass to infinity if we wish to generate 
a continuum of eigenfunctions. We take the set of eigenfunctions with positive eigenvalues as complete. 


Y^FJ^ir) (Fr(r'))*=<5-J(r-r0, 


( 20 ) 


at least for the problems of interest. For each we choose a positive Ha and take it as the frequency of the 
eigenfunction. Some additional properties of the mode functions that are required for reduction of the acoustic 
Hamiltonian to canonical harmonic oscillator form are developed in section |S.H1| of the supplementary material. 

Following similar arguments to those used in quantization of the vacuum electromagnetic field [15] (or the electro¬ 
magnetic field in nondispersive media |39j ). one can show using (20) that introducing new operators 6 a and 6]^ with 
commutation relations 


[6a, 6a'] = 0, 

&A, 6a/ = (5aa' , 


( 21 ) 


we can preserve the commutation relations (13), by expanding the fields as 

h 

A ^ ^ 


iij) = E 

A 

^(r) = 


2Ha 


6AFA(r) -kh.c., 


6.Ha 


6AFA(r) -k h.c., 


(22a) 

(22b) 


where h.c. is the Hermitian conjugate. Then, through a series of manipulations using properties of the FA(r), the 
Hamiltonian can be reduced to the canonical harmonic oscillator form 


hA = ^6Ha6],6a, 


(23) 
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where we have dropped the zero point energy which has no dynamical effect. Derivations of Eqs. (22) and (23) are 
provided in section S.IV of the supplementary material. 

Finally, returning to the physical modes by introducing 


UA(r) = 


FaW, 


^A\/p(iO 

nA(r) = -iv'p(r)FA(r) 


we can write the full displacement and momentum field operators as 


u(r) = ^^^aUaW + h.c., 

A 

7r(r) = ^^^aHaW + h.c., 

A 

with prefactors in the expansions reminiscent of simple harmonic oscillator physics. 


(24a) 

(24b) 


(25a) 

(25b) 


C. Waveguide acoustics 

The results to this point apply to any acoustic structure. We now specialize to waveguides running along the z 
direction, and choose a box of length L in that direction, which includes all x and y. Focusing on acoustic modes 
confined to the waveguide, we label the modes by a wavenumber q = 27rn/L, for integer n, and a band a which 
identifies the transverse spatial mode structure in the xy plane. Then the eigenfunctions of the previous section can 
be written 


FQq(r) — 


fa9(a:,2/)e*«"^ 

Vl ■ 


with the normalization (19) guaranteed by requiring 

J dxdyf*g{x,y) ■ taq{x,y) = 1, 


(26) 


(27) 


where the integration is over the whole x-y plane. The derived mode amplitudes (24) are then of the form 

u„,(a;,?/)e*«^ 


with 


U„,(r) = 

naq(r) = 

y) — 


Vl ’ 
Vl ■ 


1 


^ iaq{x,y), 


Llaq\fpWy) 

^aq{x,y) = -iVp{x,y)iaq{x,y). 


Note that the normalization condition (27) can be written, using (28a), as 


nlq J dxdy p{x, y) <^(x, y) ■ Uaq{x, y) = 1. 


(28a) 

(28b) 

(29) 


If we now let L —>■ oo, moving to a continuous distribution of modes, the commutation relations, Hamiltonian and 
field operator expansions respectively become 

{^aq: ^a'q^] d; 


h 

Uaq, ^a'q' 


= Saa'S{q - q'), 
























(30) 




u(r) = ^ j -^^^^^baq{t)uag{x,y) e*«^+h.c., 


(31) 


■(!')= XI / -^y-^baq{t)TTaq{x,y)e^'^^+h.C., 


while (27)-(29) remain unchanged. Here the q integrals are over the range of q for which each mode a exists, taking 
account of any modal cutoffs. Normally we work in the Heisenberg regime so the b^q are time-dependent. Note our 
convention that while u and tt are operators, the mode functions Uaq and TZaq, which carry modal index subscripts, 
are c-number quantities. 


D. Envelope functions 


To make the connection to the more familiar waveguide representation of slowly-varying envelopes, we now introduce 
envelope functions for each type of acoustic mode. We assume the excitation is centered at some wavenumber and 
factor out that dependence to produce a function varying slowly in space. However we retain the full time-dependence 
in the operators baq{t) writing 




i{q-q^)z 


(32) 


If the integral in (32) is taken to range over all q it is easily checked that we obtain the canonical equal-time continuous 
commutators 




6aa'S(z - z'). 


(33) 


In reality, the range of integration in (321 is restricted by modal cutoffs, which will temper the Dirac delta function 
in (33). However, assuming spectrally narrow envelopes far from any cut-offs, Eq. (33) should normally be an excellent 
approximation. From (31) we see that if flaq, 'Haqix,y), and •Kaq{x,y) vary little over the range of significant q we 
can approximate the field expansions in terms of envelope functions centered at q = q„: 


u(r) 


7r(r) 


cxQo J V TT 

^gzgozyM^Xu„_^^(a::,y) t)+h.c., 

CKQo 

f -k h.c. 

ciiQo ^ 

X ( 2 ^, y) <^a(2, t) -k h.C.. 

OiQo 


Corrections can be included by expanding the prefactors ^JMTa^J2^\laq{x^y) and ^h£laql‘^'^cxq{x^ y) about q = qo, 


and combining the resulting powers of {q — qo) with the exjp(i{q — qo)z) in the integrals over q to yield expressions 
involving derivatives of (l)a{z,t), but we neglect them here. 

Returning to Eq. (32), we can derive approximate equations of motion for the (f>a{z,t) in the Heisenberg picture by 































9 


expanding the dispersion relation VLaq of mode a about 


d(t)a 

dt 




ih J 


—I 




+ Vaq^ (9 - 9o) + 2 ^ago (« “ 9o)^ + ■ ■ 




„'r. j, 9^a{z,t) ^ 1.^, d^(l)a{z,t) 

^^Qgo<Pa(^,0 '^aqa aq^ Q^2 


(34) 


where 


'^<yqo 


/ 





(35a) 

(35b) 


etc. An expression for the group velocity Vaq 
supplementary material. 


in terms of the modal fields is worked out in section S.V of the 


E. Acoustic powers 


Finally, we establish some expressions concerning the power carried by the acoustic modes in the envelope repre¬ 
sentation. Classically, the acoustic power density at a point in the medium in a direction n is given by 


cl ■ ■ 


5w®(r) 

dt 


c*jfci(j.)5“(r)n' 


kl{ 


(36) 


which has the natural interpretation of power being the dot product of an applied force and the velocity of the point 
of application [35]. Using the power carried by a waveguide mode in the z direction is thus 


PtAz) = - / dxAy 


*(r) 


p{x,y) 


d^^\x,y)S^\v). 


(37) 


To construct a quantum operator corresponding to the power density w e use the symmetrized form of the non¬ 
commuting operators 7r(r) and S^^{r), in the usual way, as shown in section S.VI of the supplementary material. This 
leads to the result that the operator for the power carried by the acoustic field is 


P^{z) = E 

cx,a' ' 




(38) 


Here the contribution from the pair of modes a and o' at wavenumbers g' and q is given by 
P^a(9',9) = 


fc / O /* ^izkz 

IL / i Ln>a' 


^a'q' 


J {fl'q'{x,y))* f^q{x,y) 


h / n 


«g 


2 \ O / / 

ifl, I (jt 

2 V n^n 


(faqix^y)) 


ih / ri. 
fin 


otq 


dxAy 

dxAy 


P{x,y) 

{fa'q-i.x,y)y 

'Jpi.x^y) 


■■^Hx v)(—( 


( J 9 (fa'qf^^y) 


where we have used (281 to introduce the modified mode functions ia,q- 
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For a phonon field involving only one mode a and assuming we can neglect the q dependence of Paaiq, q') over the 
pulse spectrum we obtain the slowly-varying power operator as 



dgdg' t 

27r “9' 






(39) 


power carried by the acoustic envelope is 


Finally, it may be shown (see section S.VI of the supplementary material), that Paai^o,qa) = so that the 


P^{z) = nn^q^ Vaq, (j)\[z)(l)a{z), (40) 

and it follows that in this limit, (l)\^{z)(j)a{z) has the natural interpretation of a phonon number density operator. 


III. QUANTIZATION OF THE ELECTROMAGNETIC FIELDS 

To construct the full opto-acoustic Hamiltonian we will need similar results for the quantization of the electromag¬ 
netic field in integrated structures. The procedure is well known and we simply summarize some essential results [39) . 


A. Hamiltonian and modes 

As the fundamental quantum fields we take the electric displacement D (r) and magnetic B (r) fields with commu¬ 
tation relations 


[D\v),D^{v')\ = [B\v),B^{v')\ =0, 

[D\y),B^{y')\ = -^^5{y - y'). 

This choice, which dates back to Born and Infeld m has the advantage that the transversality of the two fields is 
easily imposed. The quantization procedure has been discussed at length elsewhere [39], and we simply quote the 
necessary. The parallels to the acoustic problem in the previous section are very apparent. 

The electromagnetic Hamiltonian operator is taken as 

^EM ^ 1 ^ b\y)B^y) dr + ^ / D\y)/3,,{{y)D^y) dr, (41) 

2Ato J 2eo J 


where 


/3ref(r) = 

eref(r) 

describes the “background” dielectric response of the waveguide structure in terms of the relative dielectric constant 
eref(r), without any acoustic effects. Note that it is straightforward to include a tensor response in eref(r) but to reduce 
cluttering the tensor notation, here we treat the material as optically isotropic. In principle, we could also extend the 
treatment to include dispersion of the dielectric [51] at the expense of considerably more complexity. For Brillouin 
processes, the linewidths of the interacting optical waves are usually narrow, and we can safely neglect dispersion 
within each optical field. 

As with the acoustic case, we are interested in optical waveguide modes with translational invariance along z. We 
find these modes {BA(r),WA} by solution of the vector Helmholtz equation 

Vx [^,ef(:E,j/)VxBA] =‘^B, (43) 

together with Ampere’s law 

DA(r) = X Ba, (44) 
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and then introduce waveguide mode functions d^k defined by 


_ d^kjx, y) ikz 

where 7 indexes the transverse spatial bands of the waveguide. We choose the normalization 

— / dxdy Prei{r) d*i,{x,y) ■ d^kix,y) = 1 , 

eo J 

and expand the displacement field operator as 

D(r)=^ J +h.c.. 

The mode operators a^k satisfy the standard commutation relations 

[cLjk ; — 0 ; 




= Sj-y'S{k — k'), 


and neglecting the vacuum energy, the electromagnetic Hamiltonian reduces to 

jjEM ^ f dkhuj^kal^k^-fk- 


(45) 


(46) 


(47) 


(48) 


B. Envelope operators 


As with the acoustic modes we can introduce envelope function operators associated with a mode 7 and a range of 
wavenumbers in the neighborhood of some kj: 




i(k-kj)z 


(49) 


The integration is to be taken over the range of wavenumbers k that we wish to associate with the center value kj. 
This allows the introduction of distinct envelope operators for fields that occupy the same spatial mode but occupy 
distinct frequency ranges (such as a pump and Stokes wave in the same mode). 

If the integrals in (49) were to extend over all k we would obtain 


,t) = 5^^>5jy5{z- z'). 


(50) 


In principle, the existence of cutoffs and the possible partitioning of each channel 7 into separate bands for pump 
and Stokes waves means that the k integrals have restricted range, but as with the acoustic fields, we assume the 
excitations are sufficiently narrow band and away from cutoff that the integrals leading to the Dirac delta function 
in (501 can be safely extended to infinity. Assuming in fact that only values of k close to kj are important for each 


mode 7 j, as we now label them, we can write 

D(r)=^ 




dk I huijk 


ajkit)d.yk{x,y)e" +h.c. 


(51) 


= E< 


/ ■^\l^^a^k{t)d^k{x,y)e^^'^ + h.c. 

■ 7 k^ix,y) J a.yfc(t)e*('""'=^)^+h.c. 


' hui^ ^ 


iJ ' 
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where we have put = oj-ykj and neglected the variation in and the k dependence of djkix,y)] as in our 

treatment of acoustic fields, corrections to these expressions can be easily identified. 

Again in analogy with the treatment of the acoustic fields, an operator for the slowly-varying part of the power 
in the waveguide can be constructed from the Poynting vector (section S.VII of the supplementary material) which 
takes the form 


7.7 ''' 

For ranges of k and k' close enough to kj, and assuming that for different 7 the corresponding kj ranges are distinct, 
we can write this as 


and since we find that 


P77 (kj, kj) = 

where vi^ is the group velocity of electromagnetic mode type 7 centered at kj, we have 

7d 


(54) 


(55) 


so that il^Az,t)'ip^j{z,t) behaves as a photon number density operator; note again the similarity with the treatment 


of the acoustic field, for which (401 is the corresponding result. The form in Eqs. (54) and (55) does not seem to have 
been presented earlier, and is derived in section [S . VII| of the supplementary materia 


Finally, in similar fashion to Eqs. (34), we find that the envelope operator obeys the dynamical equation 


d'ip'fkj (z) 
dt 


= -iui^ii^j{z) - 


dz 


+ 


^7JV 

az2 


(56) 


where 


vi^ = 


= 


dw. 


7^ 


dfc 

dfc^ 


k-k-i 


k—k.i 


are respectively the group velocity and group velocity dispersion of mode 'yj at its reference wavenumber kj. For the 
narrow bandwidths involved in SBS physics, higher dispersive terms are unlikely to be needed. 


IV. THE COMPLETE OPTO-ACOUSTIC HAMILTONIAN 


At last, we can now assemble the complete opto-acoustic Hamiltonian 


H = 


7r*(r)7r*(r) 

2p(r) 


dr+- S^^{r)A^^Ar)S'‘Ar)dr+ — 
2 J 2/io 


B\v)B\r) dr+ — 


D\r)l3^^{r)D^{r) dr. 


Being composed of different classes of oscillators, the electromagnetic and acoustic fields commute with each other. 
All quantities are taken as position dependent, varying continuously (if rapidly) across any material boundaries. Only 
at the end of various calculations we will allow them to acquire step-wise discontinuities. 

The opto-acoustic coupling is captured by the new quantity /3*'^ (r). This is the total inverse (relative) dielectric 
tensor, 


/3*7r) = (5*-^/3ref(r) -k[u (r)]), 


(57) 


which includes both the purely electromagnetic properties (the background waveguide structure) in /3ref, and the 
photoelastic and radiation pressure couplings in the correction _ Naturally, we have (r)e'^^(r) = <5*^, where 
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e^^{r) is the complete relative dielectric tensor. The coupling between sound and light enters because we assume that 
(r) depends on the displacement field u (r), through both its dependence on the strain in the material and the 
motion of the interfaces. Thus we can write 


H = H^ + + V, 


(58) 


where the opto-acoustic coupling is 


V = 


2eo 


(r; [u (r)])Zl'^ (r) dr. 


(59) 


In principle, an additional coupling arises from the dependence of the mechanical density and stiffness on the electro¬ 
magnetic field variables. These effects lead to terms quadratic in u(r) or 7r(r) (corresponding to two-phonon-single- 
photon interactions) which are of higher order than we consider here. They are also of much lower energy and the 
processes are very unlikely to be phase-matched, so they are safely neglected. 


To proceed to a set of coupled mode equations, we seek to expand the interaction Hamiltonian (59) in terms of the 


mode operators constructed earlier. The physics of the opto-acoustic interaction is introduced by writing 

/3*^(r; [u(r)]) = + S^^/3re{(r - u(r)) -f ... 


(60) 


keeping only linear terms in the strain in keeping with our neglect of two-photon interactions. The photoelastic 
tensor y) accounts for the conventional electrostrictive/photoelastic contribution to SBS. The effect of moving 

boundaries and radiation pressure enters through the second term’s dependence on the displacement u(r). Note 
that while we focus below on the effects associated with material discontinuities, this expression also accounts for 
a bulk contribution to the radiation pressure in graded index materials for which /3ref(r) varies smoothly in space. 
Equation (60) is a key expression because the symmetric relationship of radiation pressure and moving boundary 


effects follow directly from its form. This identification is what allows us to avoid the rather lengthy thermodynamic 
arguments which were required in the previous rigorous derivation of the classical coupled mode equations m- 
Expanding /3ref(r — u(r)) for small displacements we take 


/ 3 ref(i' - u(r)) « , 5 ref(r) - M (r) 


dPrelir) 


dri 


and using (57) we have for the opto-acoustic correction 

[u(r)]) =p*^'”^(x,y)5''"(r) - 


Then interaction (59) becomes 


1 


E = -^ / D\r)DJ (r) p^^^^{x,y)S‘"^{r) - . 

2eo J V \ 


dri 
( dj3re{{x, y) 


dri 


M^(r)^ dr. 


(61) 


(62) 


(63) 


Now since the phonon energy is much smaller than the photon energy the only significant terms in V will involve the 
creation and annihilation of photons. On substituting (47) for D, we normal order the photon mode operators that 


arise and neglect the resulting constant terms corresponding to vacuum fluctuation corrections to both the photoelastic 
tensor and the displacement-induced change in the dielectric properties. We thus find 




7,7' 

From and the second of © we can write 

5'-(r) = E / 5i™(r) + h.c.. 


V 


(r) dr. 


(64) 

(65) 

( 66 ) 


where S^(r) is of the form 


r,lm 


SL^(r) = 


(x,y)e 


,igz 




(67) 
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After some manipulation (see section S.VIII of the supplementary material), the interaction can be reduced to the 
form 


E 


dkdk'Aq 
(27r) 


3/2 


il/k^yk'baq j T{^k]ik'-,aq)e^^'" dz 


7,7' ,q; 

where the coupling parameter is 


+ E / J ^*hk--f'k';aq)e dz, 


r( 7 A:; ^k'; aq) = 


1 hoj^k klOrylk' Kl, 


J dxdy {d\k{x,y))* di^,^,{x,y) (^p^^^'^{x,y)s^^{x,y) - ulg{x,y) 


( 68 ) 


Note that if there is a slow variation of the nonlinear properties, due to longitudinal variation in, say, the composition 
or waveguide dimensions, then T(p/k;yk']aq) will acquire this variation too, and the integration over z in (68) would 
capture this effect. 

For an infinite homogeneous waveguide we can do the remaining integral over all z in (681 to obtain a delta function 
S{k' — k + q). Using this to eliminate the q integral in ([6^, the total Hamiltonian (58) becomes 


H — ^ ( / dk hW'yk A ^ ( / dqh^Q^q b^f^^b/yq 

7 a ^ 

+ E / a\kai'k'ba{k-k') r(7fc; ik'\a{k - k')) 

7 , 7 ' ,Q; '' 

+ E / ‘^^^^Uk-k')°^i'k'a^kT*{-ik]ik']a{k-k')). 


Observe that the final two terms explicitly display momentum conservation, and are clearly identified as describing 
anti-Stokes and Stokes processes respectively. This infinite structure form is a good starting point for investigating 
the enhancement and suppression of Brillouin scattering by adjusting the matrix elements or the density of states of 
optical or acoustic modes [52]; a simple Fermi’s Golden Rule calculation reveals much of the underlying physics, as we 
will show in a subsequent contribution |53j . Technically, the finite phonon lifetime requires the 5(k' — k q) function 
to be broadened into a linewidth function. In practice the 5 function is a reasonable approximation, with the impact 
of loss entering through the linear properties as discussed later. 


V. QUANTUM COUPLED MODE EQUATIONS 


To derive coupled mode equations for the envelope function operators, we assume that the process of interest 
destroys pump photons with transverse mode and center wavenumber (yp, fcp), creating Stokes photons with ( 7 s, ks), 
and phonons with mode-wavenumber pair {a,q) = {aa,qa)- Additional processes could be included to describe 
cascaded SBS phenomena [nisi]. We assume that the center wavenumbers and frequencies involved satisfy energy 
conservation and phase matching; that is. 


UJp — Wg -|- Qo 

kp = ks + go, 


(69) 


where we note the wavevectors can be positive or negative, and we have defined cOp = uipkp, Wg = Wgfcg, and Qo = ^^ao 9 o- 
Figure [^indicates the significance of these quantities for the processes of backward, forward intermodal and backward 
intermodal SBS. 


Then we can write (68) as 


V = dz 


dkdk'dq 


{2tt) 


2,12 pfc 




+ 1'^^ I F^ypA:; 7sfc'; ^og). 
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(xi, VL 



FIG. 2: Phase-matching diagrams for backward SBS (top), forward intermodal (bottom left) and forward 
intramodal SBS (bottom right). Black lines represent optical dispersion relations, dashed blue lines represent 

acoustic dispersion relations. 


Given the narrow bandwidths associated with SBS, we treat the coupling strengths r(7pfc; \ a^q) as constant over 
the range of wavenumbers integration: 


r(7pfc; Isk'; a^q) « r(7pfcp; 7sfcs; = F. 


(70) 


Pulling this term out, moving to the Heisenberg picture and using the definitions of the envelope functions in Eqs. (32) 
and (491, we reach 


V = r Jdz + T* J dz (l)^{z,t)tpl{z,t)tpp{z,t), 


(71) 


where tps, V’s; (!> the pump and Stokes photon fields and phonon fields respectively. Using (34) and (56), we find 
dipp{z,t) 


dt 


= ^ ^ [V'p(z,f),U] 


ih 


, , , d'ihpiz,t) 1 . d'^'ibp(z,t) r , , , ,, , N 

= -iujpipp{z,t) - Vp - ——- + -iVp - —- + —ips{z,t)(l){z,t), (72a) 


and similarly 


dips{z,t) 

dt 

d(t>{z,t) 


-iu}stps{z,t) - Vs 
—ifl^4’{z, t) — V, 


dips{z,t) 

dz 

d(j){z,t) 


1 . , d'^il)s{z,t) 

1 , ,d‘^(j>{z,t) 

a . + 


+ ■T^(l)Hz,t)'lpp{z,t), 


—i^l{z,t)^p{z,t). 


(72b) 


dt 


dz 


(72c) 
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To remove the fast time-dependence, we define 


4'p(z,t) = ^/>p(z,t)7/fia;p|up|e“P*, 

(73a) 

4's(z,t) = iljsiz,t)y'hujs\vs\e''^^\ 

(73b) 

$(z,f) = (j){z,t)^/hno\v,,\e"^‘°\ 

(73c) 


so that 

= hWp |Wp| V'p(z,t)V'p(z,t), 

= hujs ks| i^l{z,t)i’s{z,t), 

= hfl^ |li„| (j)^ {z,t)(j>{z,t). 

Thus 'hp(z, t)tE'p(z, t), 'hs( 0 ,t)'I>s(z,t), and {z,t)<^{z,t) can be identified as the power flowing in the waveguide in 
each mode at position z] this is a positive quantity with the direction along 2 given by the sign of the velocity. 
Defining a reduced coupling constant 


r( 7 A:; ^'k'] aq) = 


1 


and evaluating it at the center frequencies of our equations: 

1 


:r(7fc;7'fc';Q;(?), 


Tn = 


{Hujp) (fiuJs) IripUsU, 


=r(7pfcp;7sfcs;ao9o), 


we obtain the Heisenberg evolution equations for the field operators 


i9'l'p(z,t) i9'I'p(z,t) 1. ,9^'hp(z,t) 

dt ¥z 2 *^^ 

d^^{z,t) d'is{z,t) I. ,d^'^s{z,t) 

d^{z,t) d^{z,t) 1 . 

dt 2 *""° dz^ 


= -iujp |up|fo ^'s(z,t)$(z,f), 
= -iuJs ks|fo <^\z,t)'^p{z,t), 
= \Va \fS t)'^piz, t). 


(74) 

(75) 

(76a) 

(76b) 

(76c) 


We have thus established a fully quantum form of the opto-acoustic coupled evolution equations and shown that 
the correspondence of coupling constants for photoelasticity plus moving boundaries and the reverse processes of 
electrostriction plus radiation pressure emerge naturally from the starting point of (60). We next evaluate the 
coupling constants (or “matrix elements”) to connect them to familiar forms in the literature. 


A. The coupling constants 


We separate the coupling constants T{'yk;j'k';aq) in Eq. (74) into the contributions 

f (yfc; j'k'] aq) = fbuik(7^; I'k'; aq) + T^uTfilk; 7 'fc'; aq), 


(77) 


where 


rbuik(7^;7'^';a9) = 


eo 2^21 


V-yk'^^' VQ,q I 


dxdy y))* {x, y)p^^^'^{x, y)s^^{x, y), 


(78) 


and 


rsurf(7fc;7'^';a'3') = - 


1 


Co 2'\/2 \ Vryf^Vry/V q 


dxdy {d\k{x,y))* d\,y{x,y) Kq{x,y). (79) 


The label “surface” for Eq. (79) is perhaps too restrictive since the derivative dPref/dr^ is non-zero in graded index 
materials, however this contribution is typically weak, and in current experiments it is the surface contribution that is 
of most interest. We seek expressions in terms of the optical modes (1..^^ and acoustic modes that can be obtained 
from numerical solvers. 
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1. The bulk coupling constant 


The bulk term in (78) describes photoelasticity (in (76a) and ( |76b )) and electrostriction (in (76c)). Its evaluation 
is straightforward. Even in the limit of a step discontinuity in material parameters across a waveguide interface, the 
djk{x,y) will suffer at most a step discontinuity, as will sl^{x,y) and, by assumption, p^^^'^{x,y). Hence (78) remains 
well-defined for step discontinuities. 

Using the definition of s^^{x, y), one can show that for the standard interaction with modes (p, s, o), the bulk element 
reduces to the form 

Tbulk = 


dxdy (Ps{x, y){di{x, y))*p^^’-"^{x, y) 


-^Um{x,y) ) +iqSmz{ui{x,y)y 


2^/'^y/\VpVsV^\ Co 

which is consistent with the results in Wolff et al. m for the photoelastic coupling, there termed 


(80) 


2. The surface matrix element 


The surface term in (79) describes driving of the optical helds by moving boundaries (in (76a) and (76b)) and 
radiation pressure (in (@). Its evaluation is more subtle, since in the limit of a step discontinuity in material 
parameters j3ref{x, y) will change in a step-like fashion, but its derivative can be Dirac delta-function-like. Conse¬ 
quently, a smoothing operation is required to make sense of this term. Following Johnson et al. |55j one can show 
(see section S.IX of the supplementary material) that for a boundary contour Rc('S) parameterized by arc length s 
that separates two materials with dielectric constants e+ and e_ the surface contribution reduces to 

1 


rsurf(7^;7'^';a'?) = 


2^/'^^y\v^kVYk'Vaq\ 

X — - -) [ ['^4(R-c(s))]*dyfc-(Rc(s)) [n(s) • Uag(Rc(s))] dRc{s) 

_eo e_ e+ J 

-feo (e+ - e_) f [e||^^(Rc(s))j • (Rc(s)) [n(s) • u„,(Rc(s))] dR^s) 


■ ( 81 ) 


Here the unit normal n(s) points in the direction from e_ to e+. The full expression for rsurf(7fc; 7^^^ o:q) then involves 
a sum over all such curves separating distinct dielectrics. Note there is no ambiguity in evaluating these terms, since 
is continuous across a step discontinuity in P^e{{x,y), as is eij,^.(r). 

For the SBS combination of modes (p, s, o), the surface term can also be written as 


pMB _ 




VpVsVa 


dRc{s) (u* • n) 


1 ^ 1 


Co ve- 


1 

e+ 


(dg)*dp -h eo(e_ - e+)(n x eg)* • (n x ep) 


(82) 


which coincides with the expression of Wolff et al. m- 

A normalized form of these expressions convenient for working with numerical mode solvers is provided in section [S.X| 
of the supplementary material. 


B. Recovery of classical coupled mode equations 


Finally, the standard classical coupled mode equations 133 can be recovered from ( |76[ ) by dropping the dispersion 
terms and taking mean values for the operators: 


^(V'p) , 9(V'p) . I If ,, WM 

■Vp —— = -lojp |up| To {tps}{<P), 


dt 

di^s) 


d{4>) 

dt 


dt 
“ dz 


+ Vs 


dz 

a(V’s) 

dz 


= -iujs luslTg {(j))*{yp), 


+ v^a{4>) = |Uo| Fp (V’s)*(V'p)- 


(83a) 

(83b) 

(83c) 
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We have introduced the acoustic loss a by hand following the expression of Wolff et al. 

a=^J dxdyu*dj{riijkidkUi), 


(84) 


where rjijki is the viscosity tensor. 


VI. DISCUSSION 


Equations (76) represent a full quantum description of guided-wave optoacoustic interactions. Although it requires 


some preliminary derivations to identify the effective fields that are involved, our approach provides a direct derivation 
of both the photoelastic/electrostrictive and moving-boundaries/radiation pressure components of the SBS interaction. 
The latter has clear contributions from both surface effects at material boundaries and bulk effects due to smooth 
variation in dielectric properties. By avoiding any discussion of forces and stress tensors, the ambiguities and challenges 
of prior treatments do not arise. 

The equations of motion have a number of potential applications. In the classical limit they provide a rigorous 
confirmation of the earlier treatment m- In the quantum regime, we have a theory of opto-acoustic interactions 
that faithfully represents the photon and number statistics including any non-classical behavior. To date, quantum 
acoustic effects in guided wave systems have not been observed due to the overwhelming thermal contribution to 
the phonon field, though this may change in the near future. However, we can certainly envisage mixed systems in 
which a classical coherent state phonon held interacts with non-classical photon states in order to transfer quantum 
information between different channels, and our treatment is ideal for studies at this quantum-classical boundary. 

An obvious and signihcant extension to our work would be a complete treatment of the acoustic dissipation, which 
we plan to present in the future. To incorporate dissipation into a Hamiltonian picture, a natural approach would be 
to introduce a thermal bath of phonon oscillators rather than just the Brillouin excited mode. These modes would 
couple with the coherent modes of interest and with each other through a three-phonon collision term associated with 
anharmonicity in the phonon Hamiltonian. Tracing over the additional oscillators would lead to an effective dissipation 
on the preferred phonon mode. According to need or preference, one could derive a dissipative master equation for 
the reduced phonon density operator, or perhaps more usefully, a set of Heisenberg equations with Langevin noise 
terms associated with the loss. This kind of treatment would be particularly important in understanding the impact 
of phonon loss on the photon quantum noise. 

As mentioned earlier, another avenue is the description of enhancement and inhibition of SBS through density of 
states engineering [52] ■ In the spontaneous regime, this is well handled by a Fermi Golden Rule calculation of phonon 
generation rates [53] . as we will demonstrate in a subsequent work. 
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S.I. HEISENBERG EQUATIONS FOR ACOUSTIC FIELD 


In this section, we derive the equations of motion ([^ in section |n] of the main paper. 
Using Eqs. (50) and ([^ we can derive the relations 

1 /• [u"(r),7r*(r')] 7r*(r') ^ 1 /■ 7r*(r') [u’"(r), 7r*(r')J 


ih J 

= j 

_ 7r”(r) 
p(r) 


2p(r') 
2p(r') 


dr 

dr' + 


ih 


2p{r') 


7r*(r')(5(r-r') , 

2p(r0 


and 






ih2 J 

= --d" 
2 


(r') 


, dr , 
dr'^ ih 2 


dr'J 


dr'i 


{v'y 


dr'i 


dr'i 


dr' - L"‘ f dr' 

' dr'^ 2 J dr'^ ^ ’ dr'^ 




I d 


2 do 
d_ 
dc 
d 


= 


^ du^{v) 

dU 


1 d 


+ 2 8? 1 ■=■’"(■■) 


^ dC{v) 
drl 


.n,ki( . du^r) 

1 ) dU 


= —{c-^^\v)S^\v)), 


where in the second last line we have used ([^. 


S.II. OPERATOR AU*(r) IS HERMITIAN 


Here we show that the operator AI"^(r) of (15) is Hermitian. Consider an integral over an appropriate volume 


and assume that fields are either periodic over the volume or vanish at the surface of the volume. Then for vector 
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functions C(r) and D(r), integrating by parts twice gives 


dr 


d ({D'^{r)Y 


dr^ \ pi/2(r) 




d ( C"=(r) 


dfl 


f ( (D-(r)r 

J pi/ 2 (r) dr’- \ \dr^ \ pi/^(r) 


dr 

dr 


Now using ([^ we put = c^*"-^(r) and switching the dummy indices j O I, we can write this as 


' {D-{v)r{M-\v)C'^{v))Av = - J ^ 

_ f r (CHr))* d 


pi/2(r) dr^ 


(Minify. 


dr’ \ pi/2(r 

d f{D^(r)) 
dr’ \ pi/2(r) 


dr 


dr 


J {C’^{r)y (M'=”(r)D”(r)) dr 


and so the differential operator Al"^(r) is Hermitian. 


S.III. PROPERTIES OF MODE FUNCTIONS AND PARTNER FUNCTIONS 


Here we establish some useful properties of the mode functions FA(r) introduced in 0 to ( |^ , that are required 
to reduce the acoustic Hamiltonian to canonical harmonic oscillator form in IS.IVI 

Note that since Al"^(r) is real, if FA(r) is an eigenfunction then F^(r) is also an eigenfunction with the same 
eigenvalue wa- This may happen simply because FA(r) is purely real. In fact, as we show in S.XI it is always possible 
to choose the set of eigenfunctions {FA(r)} such that each of them is purely real. But it is often mo re con venient to 
work with complex eigenfunctions (traveling waves rather than standing waves, for example). Section S.XI establishes 
that if we include complex eigenfunctions in the set {FA(r)}, the set can be chosen so that each eigenfunction FA(r) 
is either real or, if not, there is another eigenfunction FA(r) in the set such that FA(r) = Fl(r). 

Typically the naturally chosen set of eigenfunctions will make this so; for example, if FA(r) is a traveling wave to 
the right, then FA(r) is a traveling wave to the left. We refer to FA(r) as the “partner” of FA(r). That is, each 
complex eigenfunction in the set has a partner that is also in the set. If there are purely real eigenfunctions in the 
set {FA(r)}, we take them to be their own partners. Then the set of eigenfunctions {FA(r)} is equivalent to the set 
of {FA(r)} of partner eigenfunctions, and FA(r) = Flir) for each A. Since the set of partners is equivalent to the 
original set, then from Eq. (20) of the main paper we can also write 


and 


I F*(r).FA,(r)dr = < 5 AA' 


Y^F-ir) (F“(r'))*=<5-<5(r-rO. 


(S.l) 


(S.2) 


Similarly, from Eq. (18) of the main paper we see that we have 

UA(r) = Ul(r), 
nA(r) = -fllir). 


(S.3a) 

(S.3b) 


S.IV. ACOUSTIC MODE EXPANSION 


Here we show that through use of the partner functions introduced in [S.IIH the acoustic field operators and 
Hamiltonian can be expanded in terms of the mode functions as expressed in Eqs. (22) and (23). Recall that the 

{FA(r)} are the eigenfunctions of (15) and that the |UA(r)| 


and 


{nA(r 


are defined as in Eq. (18). Both sets of 
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functions are proportional to the {FA(r)}, so we can take each to constitute a complete set of states. We can then 
expand 


^(r) = ^ 


(S.4) 


2f2A 


cfnA(r) 


where and are operators and the factors ^Jh/2U,{^ are added for later convenience. 

In the Heisenberg picture, and are time-dependent, and therefore so are u(r) and 7f(r). However, since 
u(r) and 7f(r) are Hermitian the |Ca^| are not all independent, nor are the Using Eqs. (S.3) we have 


- E (d-)' U(e ^ E (d>)' U(e 


A 

^ E "A(r). 


so that (S.41 may also be written as a sum over partner modes, 

- = E 

A 

— E 


^ cWUA(r), 


2Ha 


A"*’"*''-'- 


where we have used Oa = Ha. Then from the Hermiticity of the canonical fields we see that we require 


[e 


(i)V _^(i) 

A ^ “A ’ 




"A ; 


which may be satisfied by setting 


= &A + 4, 

cf =6 a-&^. 


(S.5) 


For real eigenfunctions FA(r) we have A = A and this just says that is (proportional to) a coordinate operator, 
and C\ is (proportional to) a momentum operator (as we will see). For partners, 6 a and 6 a are independent operators 
(or independent amplitudes in the classical case). Using (S.51 in (S.41 we have 


“W=E\/J(^A + ^y FA(r) 


(S.6) 


= E\/prW^AFA(r)+5: 


2Ha 

~h 


h 

^^a-a 




= E\E^^AFA(r)+5: 


= E 


2Ha 

AT 

W, 
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AFl(r) 


6AFA(r) -b h.c.. 
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where in the third line we use the fact that the sum is over the whole set of partner functions. Similarly 

^(r) = ^ 


(S.7) 




hil,\ 


A A 

= \/^^^AFA(r) + h.c. 




in accordance with Eq. (221. 


Postulating the commutation relations 


[A, A'] = 0, 

= <5aa' , 


(S.8) 


(S.9) 


6a, 6a' 


we find 


Xr),ff-(rO]=*^ 

A,A' 

-E 


' V Sf^A V 2 


h / 


A,A' 


' V 2f2A V 2 


h hii.A' r,t , 

- ' - b'A^bA 


bA,bl]nir) (EXl(rOr 

(FX(r)rFXl(r') 


= Y E W f E (^A (r))* Fr(r') 

A A 

= y E^aW (Fr(r'))* + y E^A W (Fr(r'))* 

A A 

= ihS^^S{r-r'), (S.IO) 

where we have used (20) and (S.2). Since [u’”(r), 7r"'(r)] = [u'”(r), 7r”*(r)] we recover the starting commutation 
relations (H). It is possible but more complicated to show that demanding the result ( |S.10[ ) one can find that the 6 a 
and b\ must satisfy (S.9). 


Now we look at the Hamiltonian in terms of the 6 a and b^. From (|16[) we have 


1 


1 


H^ = - ff*(r)f*(r)dr+ - / u*(r)X*'=(r)u'= dr. 


From the above we have 


so 


"■(■■) = E / + E ■ 

»'(r) = + < E (c;(r))‘. 

A * A * 

/'7f*(r)7f*(r)dr = - ^ EhaHa'6a6a' f FA{^)Fl,{r)dr 

A.A' 

- ^v^A^A'bi t bi, (^J nwn'W dr 

+ E ^V^A^A'bAbi, [ n(r) (n'(r))* dr 

A.A' d 

+ ^VnJhJbibA' n(r) (n,(r))* dr 
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In the last two terms, orthogonality gives A' = A. In the first two, we replace the sum over A by a sum over A and 
use the fact that FA(r) = F^(r); then orthogonality demands that A' = A. Recalling that fly = then have 

f r{r)r{r) dr = -^^;iOa6a6a- 

A A 

A A 

Then since 

M^Hr)uHr) = E \f^^lbAFUr) + E (^W)* , 


we have 


I u\r)M^'^{r)u\r)dr=J2 J n(r)n,(r)dr 

+£ (/nwn.wdr)' 

+ _5, dr 

+£viips^t(.lv(/n(r)(n^w)-dr) 

Using the same strategy as above this gives 

f u\r)M'''"{r)u'"{r) dr = ^ E ^^^aAA + h E ^^a^a^a 

A A 

+ ^ E f>-^AbAb{ + i E ^^AbibA. 


Combining (S.lI|S.lT) we have 


i?= ^E^^a (bAbi + bibA) 

A 

= E ^^A ^^aA + 2^ ■ 


S.V. THE GROUP VELOCITY 


Here we work out the group velocity of the acoustic modes in terms of the modal field providing an explicit expression 
for Eq. (35aI, a result which is needed in the next section. We take the continuous limit of Eq. (26), writing 


= VLFaqir) = Faqix,y)F'^''. 


(S.ll) 


From (17) we have 


ni^:F-(r)=M-^(r}TUr) = - 


\/p(x,y) 




d 




\Vpi^^y), 


(S.12) 


It is helpful to re-express A4”^ in terms of an operator operating on the transverse spatial variables only. 
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Applying to the mode (S.ll) gives 






5 f 


\Vpi^^y) 
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_ pnk r rk ] 

— U otq\ 1 


dr^ yfp 


(S.13) 


where the last line defines the action of the operator on faq{x,y). It follows from the Hermiticity of AI"^ that 
is Hermitian with respect to integration over the transverse plane. We can then write 

nl^f:qix,y)^Cff^^ix,y), (S.14) 


so that the f^g are eigenfunctions of £g and may be taken as orthogonal. 

Taking the inner product with and using the orthogonality of the faq, we have 

nig = j dxdy{f2gix,y))*Cfflg{x,y). 


Differentiating with respect to q gives 


on dDag 

^OlQ 1 

aq 




Since is Hermitian, we may invoke the Hellmann-Feynman theorem to simplify the right hand side: 
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Then the group velocity 
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(S.15) 


(S.16) 


(S.17) 


(S.18) 


. (S.19) 
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Swapping the dummy indices n o fc in the second term and using ([^ gives 


Vaq = I dxdy 
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= ^ / dxdy 

i l(yq 


( £n £k 

\Jaq) ^nzkz Jgq 

y/P y/P 
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fl 
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aq 


where the final line follows from (28). 


(5.20) 

(5.21) 


S.VI. ACOUSTIC POWER FLOW 


Even in the presence of coupling the displacement to the electromagnetic fields, or other forces, we expect the first 
of ([Sal still to hold. 


d , , 7r(r, t) 

= W' 

Since in general the power density at a point in the medium in a direction fi is classically given by 

du^ir) ijlm; 


Vn = -- 


dt 




the power in the waveguide in the z direction, integrated over the xy plane, is 


-Pci ( 2 :) = - dxdy 


= — dxdy 


p{x,y) 

ilifL, 

p{x,y) 


P^'’"(x,y)S''™(r) 


Azlm 


ix,y) 


du^{r) 
dr^ ■ 


(5.22) 

(5.23) 

(5.24) 

(5.25) 


where the second line follows from the symmetry properties of the stiffness tensor. 

We form the operator corresponding to the classical Pci(^) by a usual procedure. Since Pci{z) involves the product 
of the classical fields 7r*(r) and dv}{r)/dr'^^ we obtain the operator P{z) by using the symmetrized version of the 
operators corresponding to 7r®(r) and du\r)/dr"^: 


Piz) =-}- [ dxdy— . , 

^ ^ 27 ^ y(x,y) 


-'™(x,y) f , du^jr) du^jr) ^ 

' dr'^ 


dr 


^ £ ^izlm 




where we put 


^ ^ ^rplTL ^ ' 


Using (311, we see that has the form 


(S.26) 


(S.27) 
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(u^c'q' (x, ^ {yq{x, 


The first term contains parts rapidly-varying in space and time: 

and yyi{q,q') contains the slowly-varying terms, 

i^ia'iq,q) = b^a'q'baq 


-\-bai qih\^q 




X 




We write 


P[z) = pyz) + Psviz), 


where Prv{z) contains the contrib ution s from ( 9 , q') and Psv{z) those from (g, q'). Our interest is in the latter. 
Since the sums and integrals in (S.27I are over all a. a' a. and a', we n 


^ rv\-) rrom k™”, {q^ q ) ana -Tsvyz) tnose rrom q )■ uur interest is m tm 

Since the sums and integrals in (S.27) are over all a,a',q, and g', we may switch the dummy indices in the 
term on the right-hand-side of (S.29): 


(S.29) 

latter. 

second 


y^'^q’) bl,q,baq 

+baqb\,. 


Moving to normal-ordering with 


(<,,,(x,y)e*«'") {uiqix,y)e^i^)^ + (<,,,(x, y)e*«' 
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we have 


where 
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iq',q;x,y) = (<,,,(a;,y)e*«'^) «,(a:, 2 /)e*«^)^ -b «,(a;, y)e*«^) (ul,q,{x,y)e^‘>' 

T^iq) = (Kqix,y)e^‘>^y {uiq{x,y)e^^y^ + (Kq(x,y)e^^y «(x,y)e*^^)) 


(5.30) 

(5.31) 

(5.32) 

(5.33) 


The term involving T^'^{q) in Eq. (S.31) represents vacuum zero-point contributions and should give no net con¬ 
tribution to Psv{z), which is a directed quantity. Indeed, using Eqs. (28) and the property faqix, y) = 1 /))*, 

which follows from the Hermiticity of (see S.III), it can be shown that its contribution to Eq. (S.27) vanishes. 
The remaining contribution to Pgv can be written 
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Evaluating the derivatives in Eq. ( [S3^ , 

LaLiq\T,x,y) = 5^ (iq «'q')* u^g - iq' 

+ «'?')* (^^" 9 ) + (^““' 9 ') (^“ 9 ) ■ 


Using (28) we have 


^aq 


{fl^'q'T flq- ^ {ft-q'T {flq) 


^Oi' q' 


+iy/p 


y( 9 ( f^q " 

[dr‘ I 


£k 

J ol' q' 


dr'- \Vta'q'^/p 


if 

\J OlC 


Consequently, 


ptaiq^q) = 'i 


h flc 

2\ ll 


O /* ^izkz Pi O r „izkz 

ji.dy— (/:.,,)• /i + ~T 


ifi I qf 




dxdy 


ifwr 

Vp 


zkl 


A h 

Of I 


aq 




^ / /A 

ar' I Vp 


so that 


ptaiq^q) = 


Vp ^/p 


+ -g / dxdy 




ag 


yp 

- [ dxAy 
2 J 


/ A 

au i ^ 


In the second term we may exchange i and k because the other elements of the stiffness tensor are both the same to 
obtain 
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(S.35) 


which by Eq. (S.19I is simply the group velocity of the acoustic mode. The desired result (40) then follows from (39). 


S.VII. ELECTROMAGNETIC POWER FLOW 


Here we justify the relations (53) to (55) in the main paper for the optical power transport in terms of the optical 
field envelope operators. 

The operator for the power carried by the held is given by the Poynting vector which we write in the symmetrized 
form 

S(r, ^) = \ 0 ^ H(r, f) - H(r, 


Following (511, the E and H field operators are given by 


70 
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Ak ^ I hiu^k , 


E(r,t)]. 
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e*('=-*^^)^+h.c., 

(S.38) 

















































29 


where the mode functions satisfy 


^-fk{x,y) = 


^ik{x,y) 

eoe(a;, 2 /) 

1 


^'ik{x,y) = - - [V X (e^fc(a;,y)e*'"^)] 


— ikz 


ICJjkf^O 


(5.39) 

(5.40) 


It also follows from Maxwell’s equations that in lossless systems, for each mode jk, there is a partner mode 7 /c with 
k = —k, = bj^k and 


e^dx,y) = e*^{x,y) 
\k{x,y) = -h*i,{x,y). 


(5.41) 

(5.42) 


Using (S.39) and (S.40) in (S.36), the operator describing the total power flow in the waveguide is 
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The temporally slowly-varying part of this expression is 
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+ ( 6^774 } (S.43) 
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Since the k integral is over all wavenumbers including all partner modes, it is easy to show using (S.41) that the second 
term in this expression, associated with vacuum contributions, vanishes, as we would expect for a signed quantity. 

For the remaining non-vacuum contribution, since the sums and integrals are over all values we may swap the 
indices 7 ,and k, k' in the second term in square brackets to give 


pT = E 


7,7'' 


2tt 


dkdk' hujjk koj-fik 


dxdyi- e^fc x -b {e^k')* x 


n , pdik-k')-{kj-k-,)]z 




(S.45) 


where we have introduced the quantity 


p™ J dxdyz- {e^k x (h.y/fe/)* -b (eyk')* x h.yfc). (S.46) 

Finally, if the different modes 7 have very different center wavenumbers kj, then only the 7 = 7 ^ terms will contribute 
significantly to (S.45) and we may approximate 


jEM 




(S.47) 


with k) the power carried by the normalized mode functions 7 at center wavenumber k. 


S.l. Interpretation as the photon number density operator 


To convert the result in (S.47) to a simple expression involving the photon envelope operators we require the group 
velocity in terms of the fields. 


Noting that the basis functions B-yfc(r) = hjk{x,y)e'^'^^ are eigenmodes of the vector Helmholtz equation (43), the 
transverse mode functions h^k are eigenfunctions of the equation 




2 '^ik, 


where the fc-dependent operator Ok operates on a vector function f as 


Okf = V* X 
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V( X f 1- -z X z X { + ik 


z X —V( X f -b X 


and where Vt = [dx^dy,d]. It can be shown that Ok is Hermitian such that 


J dxdyr, ■ {Okh) = (^1 dxdyf 2 ■ (Ofefi) 


From Ampere’s law, we also have that 


Vt X h^^k = -i 


lfJjQUJ-y}^(l^k 


d.^ 1 . — ikz X b. 
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We now take the inner product with b*^ in (S.48) and differentiate both sides with respect to k: 

E J dxdyh*^k ■ Okhjk = E J dxdyh*^k ■ 
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where we used the normalization f dxdyh*/, ■ byk/yo = 1 which follows from (46) and Maxwell’s equations. 
By the Hermiticity of Ok, we can invoke the Hellmann-Feynman theorem to write the left hand side as 

J dxdyh*^ ■ Nfc = J dxdyh*,, • z x (z x 

+ i j dxdy b*j, • z X x hyk + b**, • V* x ^^z x 

= —2k J dxdyb*;. • z x (z x b..^^)^ 

+ i J dxdy ^b*^. ■ z x ^ - ikz x hyk)^ + (V* x b*^.) • 

= -2k J dxdy byf. • z x (z x b^/c)^ 


X b 
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+ i J dxdy {iyoU}ykd*yk + ikz x b*^.) • ^^z x byk 
= ^ J dxdy [b*j, ■ z X eyk - ■ z X b^fc] 

= z • J dxdy eyk x h*^ + e*^ x byk 


Comparing (S.52) and (S.53) yields 

dtOyk 


dk 


= • y dxdy Byk X h* J, + e*^ x byk 


HUJ-yk 


P™{k,k), 


Finally, from (S.46) we then havey™(fc, A:) = hioi^vi^, and from (S.47) with (49) we obtain (55) 

P™(0) K^hwi^V^y 'ip\^{z)'lpyj{z), 
id 

1 exact analogy with the acoustic result in (40) but allowing for the sum over electromagnetic modes. 


(S.53) 


(S.54) 


(S.55) 


S.VIII. SIMPLIFICATION OF THE OPTO-ACOUSTIC INTERACTION TERM 


Here we show how the interaction (64) may be reduced to the form shown in (68). Inserting the expansion of the 
strain tensor (66) into (64) gives 


7 ,7',a 

X J {dyk(.x,y))* di^,k>(.x,y) ^y*^'’"(x,y)s(,™(x,y) - u^aq(.x,y) ) dxdy 


1 

H— 
eo 


7,7 ,CX 


X J {d\k{x,y))* d‘y,y{x,y) ^y*^'’"(x,y) (s^™(x, y))* - {u^q{x, y))*'^ dxdy. 
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Since the inverse dielectric tensor is symmetric even under strain, we have y) = y), and swapping the 

dummy indices k, k' in the second term gives 


^ J dkdk'dqal^^ttYk'baq] 


hU^k klO^'k' kTlgg /" 

4tt \ 47r V 47r / 


eo —' J ' V 47r V 47r V 47r y 

7 , 7 ' ,Q: 

X f id’yk{x,y))* di^'k'i^^y) (-^^^0^^'^aqix,y)]dxdy 


1 hwjk 

/ hijo^'k' 

/hflaq f 

V 47r \ 

1 47r ' 

\/ dTT 7 ® 


X J (di^'k'i^^y)) d\k{x,y) (^p^^^'^{x,y) {s^^{x,y))* - 6^^ iu^q{x,y))*^ dxdy. 


We can now write this as 

E 


E J <^lk^yk'baq J r{jk;yk'-,aq)ed'^' dz 

7,7',Q: V '^) 

+ E / blga\,^,a^k J T*{-fk;-f'k';aq)e-d'^'-'^+'i^^ dz, 


(S.56) 


where the coupling is characterized by the slowly-varying coefficients 

r(7fc;yfc';ag) = 


X / dxdy ((i!yfe(x,y))*d^,fc,(x,y) y)s^™(x, y) - (5*^ f d^ref{x,y) \ , 


S.IX. SMOOTHING THE SURFACE MATRIX ELEMENT 

In the (x, y) plane we can in general identify a number of curves C that indicate where Preiix, y) will change 
discontinuously from one value to another. These may be straight or curved lines. We only contemplate discontinuous 
changes in j3^^f{x,y), adding up the neighborhoods of all such curves identifies the regions where Prefix, y) is assumed 
to vary in the xy plane. We write R = ix,y), and parameterize such a curve by Rc(s) = (xc(s),yc(s)), and for a 
given curve let s range from 0 to 1. As s increases along the curve we have 

ATi ( \ -dxc(s) .dyc(s) 

dRc(s) = X—,-ds -I- y — ,-ds 

ds ds 


where the length 


= u(s) di?c, 


dRr = dsi 


dxc(s)y /dyc(s)y 


(S.57) 


and the unit vector 


u(s) = 


We introduce a normal to the curve as 


^ drCcfs) I ^ d'Vc(s) 




(S.58) 


n(s) = z X u(s) 


^ daic(s) I ^ d-i/cC-s) 
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FIG. S.l: Geometry for smoothing of discontinuous fields at waveguide interfaces. 


and in the neighborhood of the curve we can specify the points in the xy plane by (s,C)) where 

R = Rc(s) + Cn(s), 


(S.59) 


or 


X = Xc{s) + c (x • n(s)), 

y = 2/c(s) + C(y • n(s)) • 

For fixed s, /3ref changes as ^ passes from < 0 to > 0; that is, it is only a function of We assume now that the 
change in j3reiix,y) occurs only in such a small region about R = Rc(s) (we will eventually take that change to be a 
Dirac delta function there) that the mapping from (C, s) to {x,y) is one-to-one. Then we can write 


/ dPrefjx, y) \ 

\ dr^ ) 




{x,y) 




n(s) • Uag(a;(s,C),y(s,C)), 


(S.60) 


and we have 


- J dxdy {d\^{x,y))* dl^,^,{x,y) ^ u'^q{x,y) 

= -J |J(s,C)|dsdC {d\kix,y))* dl^,f,,{x,y) n{s) ■ Ucq{x, y), (S.61) 

where in the second line we understand x = x{s,(^) and y = 2/(s, C); that is, we have switched integration variables 
from X and y to s and C. The Jacobian 


d(s,C) = 


dx dy 
ds ds 
ox oy 

ac ac 


f dxc{s) , an(s) 
■ 




\ ds 9s 

Now at each point (s, (") we can use the local reference frame to identify 

{djk{x,y))* dl^,f,,{x,y) = {d:^k{x,y))* d^,f,,{x,y) + y)) • y). 


X • n(s) ) . 


(S.62) 


(S.63) 


dYk'ix,y) = n(s) • dYk'{x,y), 

dyfe/ (x, y) = 22 .- dy fc/ (a;, y) -k u (s) u (s) • dyk' (x, y) 
= dyfc'(x,j/) - n(s)n(s) • dYk'{x,y), 


where 
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again recalling that x = x{s,Q and y = y{s,Q- Then we can write (S.61) as 

- J dsdy {d\f,{x,y))* dyk>{x,y) ^ u'^q{x,y) 

= - j |T(s,C)|dsdC {d^^{x,y))* d^,f,,{x,y) «(«) ’ ^aq{x,y) 

- J\Jis,C)\^sdc(d^\f,{x,y)j ri(s) •u„q(x,y). (S.64) 


From the relation 


follows 


^ref(C) “ 

d/3ref(C) 


1 


PreliC) 

1 deref(C) 


dC e?ef(C) dC 

Now the simplest characterization of the variation of /3ref(C) would be to write 

p,ei{C) = P- + {P+-P-)d{0, 


(S.65) 


(S. 66 ) 


(S.67) 


where 9{() is th e step function, /3_ is the value of /3ref(C) for negative and /3+ is the value of /3ref(C) for positive (. 
Similarly, from (S.65) we can write 


1 


1 1 


£ref(C) - ^1 ^(0- 


(S. 68 ) 


To differentiate with respect to ( and then integrate in (S.641 we smooth these functions. We introduce a smoothing 
function gi{C) which is non-negative, peaked at ^ = 0 , satisfies 

y’ 5 /(C)dC=l, (S.69) 

for each I, and approaches a Dirac delta function as Z —>■ 0. Then for finite I we have smoothed functions 

^ref(C)= J 9l{C-C)PrAC)dC. 
eref(C) = j 9l{C - OSrefiC) dC'- 

One strategy for evaluating dp^efiO/dC is to take 

d/3ref(C) , d^ref(C) 


dC 


Alternately, using (S. 66 ), we could take 


d/?ref(C) 

dC 


=iP+-P-)9i{C)- 

1 d£ref(C) 

£?ef(C) dC 


(S.70) 


(S.71) 


4f(C) \P+ P- 


Using ( |S.70[ ) and ( |S.71[ ) in the two right-hand expressions of (S.641 respectively, gives 

uL{x,y) 


dxdy {dl^^{x,y))* d\,^,{x,y) 


= -{P+-P-) J |J(s,0|dsdC {d^k(.x,y))* d^,^^,{x,y) 9 i{ 0 n{s) ■Uaq{x,y) 


+ ( ^ - ^ ) / l•/(s,C)|dsdC■ 


4f(C) 


gi{C)n{s) ■Uaqix,y), (S.72) 
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where we still understand x = x(s, (} and y = y{s, O- Now we can let Z —> 0 in both terms, because the rest of the 
integrands are continuous about C = 0- Recalling (S.69), we have 


( d(3rei{x,y) \ ^ fc 


- J dxdy {dl^^{x,y)) d\,^,{x,y) J ^ag(x,y) 

-(/3+ - I3-) J |./(s,0)| {d^,,(Rc{s))y d^,,,,{'Rc{s)) n(s) • u„ 5 (Rc(s))ds 

+ y I J(s, 0)1 (e1^,(Rc(s))) • ei|,^,(Rc(s)) n(s) • u„,(Rc(s)) ds. (S.73) 

where we have used the fact that x{s, () —)■ Xc{s), and y{s, C) —> yds), as C —>■ 0. Finally, we have 


\J{s, 0)1 ds = 


f dxc(s)\ - ^ f dyc(s) 1 . .. , 

' y • n(s) - —- X • n(s) 


\ ds 


ds 


ds 


f dxc{s)\ f dyds) 


V V 

the element of length along the curve. So we can write 


ds 


f dl3ref{x,y) \ ^ fc 


ds = di?c(s), 


dxdy [ddix,y)] ddkdx,y) J '^aqix,y) 

-{13+ - I3-) J [(i:J:fc(Rc(s))]*(i:^,fe,(Rc(s)) [n(s) • u„,(Rc(s))] dRds) 

+ J (e"fc(R-c(s))) ■ei|,fc,(Rc(s))[n(s)-u„g(Rc(s))] di?c(s) 

= (^ - y^) j {dykC^ds)))* dykd'^ds)) [n(s) • u„g(Rc(s))] di?c(s) 

+ 4 (^+ - £-) [ (e 7 fe(^c(s))) • ei|,j,,(Rc(s)) [fi(s) • u„g(Rc(s))] di?c(s). 


(S.74) 


from whence (80) follows. The full expression for FgurfCT^; 7^^^! Q^?) then involves a sum over all such curves where a 
transition from one dielectric constant to another occurs. Note there is no ambiguity in evaluating these terms, since 
d:|^(r) is continuous across a step discontinuity in /3i.ef{x,y), as is e^^.,(r). 


S.X. REDUCED MATRIX ELEMENTS 


Using the normalization conditions (29) and (46) we can write the matrix elements in the form 
f (yfc; aq) = 


'dlX'ykX'yfk' '^aq I 

Jdxdy (d;fc(x,y)) d^^,^,{x,y)p'^^”"{x,y)S^^{x,y) 


/dxdy/3ref(r)d*fc(x,y) ■ A^k{x,y) [j dxdy p{x,y)n%d^,y) ■\iaq{x,y)] 


1/2 


^8- 8^' 


(S.75) 


+ (e+ - £_) Y 


/( 

djkdds))) 

* 

dYk'dc(s)) (n(s) ■ Uag(rc(s))dRds) 

/ dxdy/3ref(r 

41 

)d*j,(x,y)-< 

(4kdds)) 

i7k(x,y) 

)'- 4 . 

1 /2 

[/dxdyp(x,y)u*^(x,y) • Uag(x,y)] 

rc(s))n(s) • Ua,(rc(s))di?c(s) 

f dxdy/3rei(r)d*^(x,y) ■ 

ijk (^5 y) 

1/2 

[/dxdyp(x,y)u*^(x,y) • Ua,(x,y)] 


The advantage of this form is that it can now be used regardless of how the mode fields are normalized. Again, the 
last two lines should be summed over all curves that contribute. 
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S.XI. THE ORGANIZATION OF EIGENFUNCTIONS 

This section establishes the basic properties of partner eigenfunctions for a Hermitian operator that are invoked in 
section IS. fill 

Consider a Hermitian operator, schematically H{x, the eigenvalue equation is 

H{x,-^,...)f{x) = Xf{x). (S.76) 

Hermiticity guarantees real eigenvalues and the fact that eigenfunctions of different eigenvalues are orthogonal. The 
inner product of two such functions vanishes, where the inner product of g(x) with f{x) is 

J 9 *{x)f{x)dx (S.77) 

We consider eigenfunctions that are normalized, so 

I r{x)f{x)dx = l. (S.78) 


We want to consider first a number of degenerate eigenfunctions, all with the same eigenvalue. Suppose now that 
besides being Hermitian, H is also real. Then if f{x) is an eigenfunction, f*{x) will also be an eigenfunction with the 
same eigenvalue. 


H{x,^,...)r{x) = xr{x) (S.79) 

For a given f{x), of course one possibility is that f*{x) is just a constant phase factor times f{x). Then f{x) could 
be readjusted to be purely real (or purely imaginary), for example. 

Suppose this is not the case. Then f{x) and f*{x) are linearly independent, and they span a two-dimensional space. 
Of course, they need not be orthogonal. That is, there is no guarantee that the inner product of f*{x) with /(cc). 


J f{x)f{x)dx 


(S.80) 


vanishes. If it does, we call f{x) and f*{x) “partner” eigenfunctions. Suppose now that 
can construct partner eigenfunctions from f{x) and f*{x) in the following way. 

First find 


(S.80) does not vanish. We 


c{x)=N{f{x)+r{x)), (S.81) 

where A is a real normalization constant; c(x) does not vanish, because by assumption f*{x) is not just a multiple 
of f{x). If we choose N to be real, then c{x) is also purely real. Now take out from f{x) the amount proportional to 
c(x), 


f{x) = f{x) - c{x) J c[x')f{x')dx\ (S.82) 

where we do not need c*{x') in the integral because c{x') is real. Of course f{x) cannot vanish everywhere because 
otherwise f{x) would just be proportional to c{x) and then f{x) would just be a phase factor times a real function. 
Now by construction c{x) is orthogonal to f{x), 


c{x)f{x)dx = 0. 


(S.83) 


Perhaps f{x) is purely real; if so, normalize it and call the result s{x). Perha ps f{x ) is purely imaginary; if so, divide 
by i, normalize it and call the result s{x). If f{x) is neither, note that from (S.83) we have 


c{x)f*{x)dx = 0, 


(S.84) 
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since c{x) is purely real. Then 


f{x)+f*{x) (S.85) 

is a real function that is orthogonal to c{x); it cannot vanish everywhere because we have assumed that /(x) is not 
purely imaginary. Now normalize this function and call it s(x). 

Whatever route we have taken to get s(x), we now have two real functions c{x) and s(x) that are orthogonal to 
each other and normalized, 


/c=(x)dx = l, 
/.^(x)dx = l. 
c{x)s{x)dx = 0. 


They span the space spanned by /(x) and /*(x). We can then form partner functions for this subspace, 

fi{x) = ^(c(x) +*s(x)), 
fiix) = ^{cix) - isix)). 

These functions are normalized, 


and orthogonal, 


J [fiix)]* fi{x)dx = 1, 

j If^x)]* f*{x)dx = l, 

J [fi{x)]* fi{x)dx = 0, 

J [fi{x)]* fi{x)dx = 0. 

So we have constructed partner eigenfunctions /i(x) and /i (x) that span the subspace spanned by /(x) and f*{x). 
Suppose now there are more eigenfunctions with the same eigenvalue, which are orthogonal to /(x) and /*(x). Call 
one of them g(x). Then g(x) must be orthogonal to /i(x) and /f(x) since they span the same subspace as /(x) and 

Fix), 


J 9 *ix)fi(x)dx = 0, 
J 9 *{x)f*(x)dx = 0. 


(S.86) 


Now if g{x) is an eigenfunction of H{x), then g*{x) is an eigenfunction of H{x) with the same eigenvalue. Suppose 
g*{x) is not just a constant phase factor times g(x); then g*{x) and g{x) span a two dimensional subspace that, since 
from (S.86[) we have immediately 


j g{x)fi{x)dx = 0, 
J 9 {x)fiix)dx = 0, 


has no overlap with the subspace spanned by /i(x) and /i (x) . So from g{x) and g*{x) we can form two partner wave 
functions / 2 (x) and /^{x) that are orthogonal to each other and each orthogonal to each of /i(x) and /i (x). 

Thus we can proceed and organize our eigenfunctions. As we investigate all the eigenfunctions of a particular 
eigenvalue we will sometimes hnd it is possible to immediately make an eigenfunction real (as we could have, for 
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example, if g*{x) had simply been proportional to g{x) with a constant phase factor), or otherwise we can establish 
partners. So we can imagine listing all our wave functions grouped in the following manner, 

h{x) fl{x) (S.87) 

h{x) /a (a:) 

h{x) fsix) 


fN{x) f^ix) 
fiix) 
fii{x) 


Here the Roman numerals indicate real wave functions that “don’t have partners”; we take them to be purely real. Of 
course, if we have an even number of real wave functions without partners we can start combining them into partners. 
For example, in the list above we could replace fi{x) and fii(x) by the partners 


fN+i{x) 

f*N+i{x) 


^ i.fi{x)+ifii{x)), 
■^Ui{x)-ifn{x)). 


If we have an even number of eigenfunctions of a particular eigenvalue, then we could pair them all up in partnerships. 
If we have an odd number then there must be at least one “unpartnered” wave function. It is also possible to “divorce” 
some partners; for example, in place of fz{x) and f^(x) we could choose the real functions 


C3(a;) = ^ ifsix) + fU^)) ^ 

S3{x) = -^ (faix) - /a (a;)). 


But it is often convenient and natural to have wave functions in partnerships. In any case, we assume that we have 
eigenf unctions organized according to (S.87). However, we henceforth write fi{x) as fi{x), and so on, so the list 
(S.87) can be given as 


fi{x) h{x) (S.88) 

f2{x) hix) 

h{x) h{x) 


fN(x) f/fix) 
fiix) 
fll{x) 


Then if we denote a general eigenfunction by fa{x), the list of possible Js is 

l,I,2,2,3,3,...iV, IV, 7, II, III... 


(S.89) 


These eigenfunctions are all orthogonal, 


j fai.x)fa'{.x)dx = 5aa'- 

as a and a' range over this list. Associated with a list of as we introduce a list of ds, 

1,1,2, 2, 3,3,... IV, N, I, II, III... 


(S.90) 


(S.9I) 
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That is, if a is one of a partnership, a is the other partner; if a is a real wave function, a is that wave function itself. 
Clearly 

(S.92) 

a cx 

and 

f*{x)=Ux), (S.93) 

either because a identifies the partner of a, or because faix) is real, in which case fa{x) can be considered its own 
partner. 

Now if we consider the eigenfunctions of a whole range of eigenvalues A we can do the same sort of organization 
within the subspace of each eigenvalue. Then we can let a range over the whole set of labels of all eigenfunctions of 
all eigenvalues. For a given a we identify the eigenvalue by Aq,. Then over this whole range of as we have 

J fai^) fa'ix)dx = Saa', (S.94) 

where between eigenfunctions associated with different eigenvalues the orthogonality holds because of Hermiticity of 
the operator, while between eigenfunctions associated with the same eigenvalue the orthogonality holds because of 
the construction we have adopted. We still have generally 

fUx)=fa{x), (S.95) 

and of course 

Aa=Aa. (S.96) 




